RATE OF CONVERGENCE IN NONLINEAR HARTREE DYNAMICS WITH 

FACTORIZED INITIAL DATA 



LI CHEN AND JI OON LEE 

Abstract. The mean field dynamics of an A-particle weekly interacting Boson system can be de- 
scribed by the nonlinear Hartree equation. In this paper, we present estimates on the 1/N rate of 
convergence of many-body Schrodinger dynamics to the one-body nonlinear Hartree dynamics with 
factorized initial data with two-body interaction potential V in L 3 (R 3 ) + L° 



1. Introduction and Main Result 

1.1. Setting and history of the problem. In a non-relativistic case, the dynamics of an N- particle 
system is governed by the Schrodinger equation. For a real physical system, however, N is usually very 
large so that it is hopeless to solve the TV-body Schrodinger equation directly. There were many efforts 
to describe such a system by approximating by a simpler dynamics. One of the most important cases 
is a system of TV-weakly interacting Bosons, which can be approximated well by using the nonlinear 
Hartree equation. 

We consider a system of iV-interacting three-dimensional Bosons in ]R 3 , described on (L 2 (R 3Ar )) s , 
the subspace of L 2 (M. 3N , dX^) consisting of all symmetric functions, where Xm '■— {xi,X2, • • • ,Xn)- 
Given the two-particle interaction V, the mean-field Hamiltonian of this system is 

N N 

H N :=-^^ j + - y £ i V(x i -x j ) ) (1.1) 

3=1 *<3 

where Xi e R 3 are the positions of the particles and Aj denotes the Laplacian with respect to i-th 
particle. Note that the coupling constant 1/N guarantees that the kinetic energy and the potential 
energy are typically of the same order. 

We assume that the initial state is factorized, i.e. 

JY 

1>(X N ) = l[<p(x j ) (1.2) 

3=1 

for some ip € iJ 1 (M 3 ) with \\(p\\ l 2 (m. 3 ) = !• Its time evolution, ip t , satisfies the A^-body Schrodinger 
equation, 

id t ipt = H N ip t (1-3) 

with the initial data ?Ao = V J - 

It is well known in this system that the factorization is approximately preserved under the time 
evolution, and, in fact, we can expect that 

N 

M*n)^ l[vt(xj), (1.4) 

3=1 

where (ft is given by the solution of the nonlinear Hartree equation, 

id t tp t = -Ap t + (V * \(p t \ 2 )<p t , <ft\t=o = <P- (1-5) 
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To give a meaning to (|1.4j) . we introduce marginal densities. The marginal density associated with 
ipt is dehned to be the orthogonal projection onto ipt, and we denote this by bracket notation, 

7t = KWt|- (1.6) 

The kernel of j t is given by 

lt (X N ; X' N ) = i> t {X N )Wt{X' N ). (1.7) 

(k) 

For k = 1, • • • , N — 1, we also define fc-particle marginal density 7 t associated with ipt by taking the 
partial trace of jt over the last (N — k) particles, i.e., the kernel of 7^ is given by 



O) 
It 



{ x k]X' k ) = J dx k+ \ ■ ■ ■ dx N j t {X k ,x k+ i, ■ ■ ■ ,x N ;X' k ,Xk+i, • • • ,x N ), (1.8) 

where X k = (xi, ■ ■ ■ , x k ) and X' k = (x[, ■ ■ ■ , x' k ). 

In terms of marginal densities, (jl.2p is equivalent to 

7 = 7o = |^)<^, (1.9) 

and it is expected to show that its time evolution satisfies the convergence 

7 t ( *> -> \<pt)(<Pt\ 9k as N -> 00 (1.10) 

in the trace norm topology. 

Spohn [IT] proved (fT7T0|) . assuming the interaction potential V is bounded, and Erdos and Yau [3] 
extended this result further to obtain a rigorous derivation of the Hartree equation for the Coulomb 
interaction case. These results were based on the study of dynamics of marginal densities, which is 
governed by BBGKY hierarchy 

^d k) = E [ - A„ 7t (fc) ] + ^ E [V(xi - x,),7 t (fe) ] (1-11) 



A^l "J' it J 1 N ^ 

j'=i *<i 



N — k 1 



E Tr fc+i - Xfc+i), 7 t ( 



j=i 

where Tr/j+i denotes the partial trace over the (fc + l)-st particle. In this method, compactness of the 
sequence of marginal densities {"ft }k=i 1S nrs t shown, then that any limit point of the sequence is a 
solution of the infinite hierarchy is proved, and finally, proving the uniqueness of the solution gives the 
desired result. 

Another approach was introduced by Hepp [5] and extended by Ginibre and Velo [H [S]. In this 
approach, the time evolution of coherent states was studied in the second quantized Fock-space repre- 
sentation. Using this method, the rate of convergence in (|1.10[) was proved by Rodnianski and Schlein 
|10j . where they proved that, for factorized initial data and for the interaction potential such that 
V 2 (x) < C(l - A x ), 

Trhf - \<Pt){<Pt\ 9k \ < %r (1-12) 



N 

for some constants C and K . 

Knowles and Pickl [9] introduced a new method to prove (|1.12[) where the interaction potential 
V E L 2 (R 3 ) + L°°(M. 3 ). They also extended the result to more singular interaction potentials, where 
the rate of convergence becomes weaker as the singularity in the interaction potential stronger. 

The optimal result for this quantitative estimate must be of the order 1/N, since the bound of the 
form 

Tr|7 t (fc) -|^)(^r|<^r (1.13) 
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is proved for the coherent case by Rodnianski and Schlein [10] and for the factorized initial state case 
with the bounded interaction potential by Erdos and Schlein [2] . 

In this paper, we first observe that the IA2 dynamics, which was introduced by Ginibre and Velo 
[4, 5 , gives the 0(1 /N) rate of convergence. The goal of this paper is to show that for factorized initial 
data, there exist constants C and K and such that 

Tr| 7t (fc) - |^ t H < — (1.14) 

for V G L 3 (R 3 ) +L°°(R 3 ). 

1.2. Notations and tools to be used. Let Fx be the Fock space of symmetric functions, i.e. 

:=0(L 2 (]R 3 ")) s , (1.15) 

17 >0 

where (L 2 (R°)) s = C. Here, s denotes the subspace of symmetric functions. A vector %j) in Fx is a 
sequence ip = {V^ }n>o of n-particle wavefunctions G (L 2 (R 3 ")) ;j . The scalar product between 
ipa,ipb G Tx is defined by 

n>0 

and we will omit the subscript Fx from now on. We let 

O:={l,0,0,---}G F x , (1.17) 
which is called the vacuum. We will also make use of the space 

F:=L°°([0,T],F x ). (1.18) 

For n > 1, let 

n 

F n := {Ut,x n ) l g Lf^xJ* n f| /.?/.'; ,/.;., /.\. ,, ., ,| (1.19) 

with 

ll^||F„ = ||en(i,^n)IL r i= +, max {||^(t,X„)|| L?i 6 L » £ " 1, (1.20) 

and let Fq — C. Here, we again consider a fixed time interval t £ [0, T] only. We define another space 

F~@F n , (1.21) 

n>0 

equipped with the norm 

II • ll| = EH "Ilk- (1-22) 

n>0 

On Fx, the creation operator aj. and the annihilation operator a x for x G R 3 are defined by 

1 " 

( a l^) (n H x i,--- ,x n ) = —=^28{x - Xj)^ 71 '^ {x 1: - ■ ■ ,Xj-i,x j+1 ,--- ,x„), (1.23) 

(a^) (n) (^ir-- ,x„) = v^+l^ (n+1) (a;,x 1 ,--- ,x„). (1.24) 
For / G L 2 (R 3 ), a^(/) and a(f) are given by 

« f (/) = / dx/(x)4. (1-25) 



a(f) = / dxf(x)a x . (1.26) 
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The creation operator a'(f) is the adjoint of the annihilation operator a(f), and they satisfy the 
canonical commutation relations 

[a(f),a\g)] = {f,g)i*0P), [«(/), a(g)] = [ot(/), J(g)} = 0. (1.27) 
The number operator AT on Tx is defined by 

M := \ dxa\.a x , (1.28) 



and it also satisfies 

(A/V>) (Tl) =W> (n) - (1-29) 
We define the Hamiltonian on Tx by (H-Nipp™' = Hnip - Using the operators a x and a\, it can 
be rewritten as 

H N = J dx a x (-A x )a x + J J dxdy V{x - y)a x a\a y a x . (1.30) 

For / 6 L 2 (R 3 ), the Weyl operator W(f) and the coherent state i/j(f) are given by 

W(/):=exp(at (/)-«(/)), (1.31) 

and 

V>(/) := W(f)Q. (1.32) 

The Weyl operator W(f) satisfies 

W(J) = e- |l/l|2/2 cxp( a t( / )) cxp (_ a (/)). (1.33) 
Important properties of Weyl operators are collected in Lemma 2.2 of |10) . 

1.3. Main theorem. In this paper, for factorized initial data and the interaction potentials V € 
L 3 (R 3 ) + L°°(R 3 ), we extend the convergence rate 0(1/Ni) in Rodnianski and Schlein [TU] to 0(1/N). 
More precisely, we prove the following theorem. 

Theorem 1.1. Suppose that V = V\ + V 2 where V\ € L 3 (R 3 ) and V 2 £ L^QR 3 ). Let 7^ be 
the one-particle marginal density associated with the time evolution of the factorized initial state 
{0, • • • , 0, ip® N , 0, • ■ • } as in (|1.8p . Then there exist constants C and K, depending only on 
||Vi||l3, and ||p2||i°° such that 

Tr 1% t -\<Pt){<Pt\ <—■ (1-34) 

Remark 1. Using the same arguments used to prove (|1.34[) , we can extend the result to the higher 
marginals 7^ t- 

Remark 2. For the proof of Theorem \l.l\ we first assume that t € [0, T] and prove the bound 

m Ce*r 



W,t-M{Vt\ <^^- (1-35) 



Tr 

Letting T — t, Theorem \1.1\ follows. 

Remark 3. What we actually prove in this paper is 

Ce Kt 

7 g -\<Pt)(<Pt\ hs <—, (1.36) 

which implies f|l .34[) . Here, \\ ■ \\hs denotes the Hilbert- Schmidt norm. The argument for this implica- 
tion can be found in Remark 1.4 of |10) . 

Remark 4. Note that the nonlinear Hartree equation (jl.5[) is globally well-posed when V satisfies the 
assumptions of Theorem ] 1. 11 See Remark 1.3 of [10] for further detail. 
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In [10], the authors had proved the convergence rate is 1/N with initial data is coherent state. Then, 
by writing the factorized state 



{0, • • • , 0, <p* N , 0, • • • } = (flt( ^- )A ' n = d N ^- e * eN W(e-* e ^ 

VN\ Jo 27r 



)n (i.37) 



with the constant 



they could prove that the rate of convergence is of order 1/y/N. 

To explain our idea more explicitly, we need first to give a brief outline of the proof of the main 
theorem in [TO] . 

1.4. Outline of the idea in [10] . Let F$Jx,y) be the kernel of the one particle marginal density 
associated with the time evolution of the coherent state W(v~N(p)Sl, which is defined by 

r$,(z,v) = ^{e- mNt W{^N^,ala x e- mNt W{^N^), (1.39) 

We expect that the limit of the kernel of one particle marginal density is Tp t (x)<pt(y), thus we expand 
r^ t (a;, y) in terms of (a x — V~N<fit(x)) and (at — \fNTpi(y)). Then, we get 

T N,t( x >y) = pt&viiv) 

+ i(fi, W\VN V y nNt (al - ^/NJpl(y))(a x - VN<p t (x))e- mNt W(VN<p)fl) 
+ f^-(n, W^{VN<p)e mNt (al - \fNTp- t (y))e- mNt W(^/Nip)n) (1.40) 



N 

+ ^}(n,W\VN<p)e mNt {a x - VN<p t {x))e- mNt W(VN<p)n). 
v N 



It was shown by Hepp [6] (also by Ginibre and Velo [1J[5]) that 

W*(VN(p a )e <H »<- t -')(a x - VN(p t (x))e- mN ^W(VN(p s ) =tf(t;s)a x U(t;s), (1.41) 
where U(t; s) is a unitary operator defined through 

U(t- s) := W\VNcp t )e- inN ^-^W(VN(p s ). (1.42) 
The operator U(t; s) also satisfies 

id t U(t;s) = (H 2 {t) + H 3 (t)+H i )U(t;s), U(s;s)=I, (1.43) 
where the generators H 2 , H 3 , and H4 are defined as follows: 

H 2 (t) := / dxV^aJ.V^a^ + / * |^i| 2 )(x)aJ.a ;E + // dxdyV (x - y)Tft(x)ip t (y)a y a x 



+ \ 1 1 dx dyV(x~y)(ip t (x)ip t (y)a x al + ipt(x)ipt(y)a x a y ) 7 (1.44) 



-ff 3 (t) := -7= / / dxdyV(x - y)tpt(y)alal<i x + —= / / dxdyV(x - y)ip t (y)a x a y a Xl (1-45) 



= J J dxdyV(x - y)vt{v)aW y a>x + 7^ 

^ 4 := 2N I dxdyV ( x ~ y) a l a l axa v ( 1 - 46 ) 
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Using U(t; s), we can rewrite (|1.40[) as 

T%{x,y) - ipt{x)Wt{y) (1-47) 

1 ^^.<n,Mt( t; o) tW(t;0)n> + ^}(n,U^{t-,0)a x U(t;Q)n). 



'N " V^V 
Proposition 4.1 of 4 shows that there exists a unique group of unitary operators U 2 (t; s) satisfying 
id t U 2 {t- s) = (H 2 (t) - #o)W 2 (i; s), W 2 (s; s) = 7, (1.48) 
where H^if) is defined by 

H 2 (t) :=e lHot H 2 (t)er tHot . (1.49) 

We define 

U 2 (t; s) := e- lHot U 2 {t; s)e lH ° s . (1.50) 
We note that U 2 is well-defined on Txi since the free evolutions e~ lHot and e lH ° s are well-defined on 
Jx, and the unitary operator U 2 (t; s) is also well-defined on Fx- Since the operators U 2l e~ lHot , and 
e jffoS are bounded in V(J\f 5 ) for any 8 E R, so is U 2 . 

Furthermore, it can be proved from Proposition 2.2 of [5] and from Lemma 1 7 . 1 1 that U 2 is strongly 
diffcrentiable from Q(H + A/"), the form domain of H + Af, to its dual Q* (H +Af). Thus, on a dense 
subset T)(Ho) (~l T>(M) of Txi we can check from the definition that 

id t U 2 (t;s)=H 2 {t)U 2 (t;s), U 2 (s;s)=I. (1.51) 
For simplicity, we will use notations 

U(t) :=W(t;0), W 2 (t):=W 2 (i;0). (1.52) 

To estimate the difference between r^ t (x, y) and ipt(x)lpi(y) as in (ll.47p . one needs to find a good 

approximation of W(i; 0), which conserves the parity of the number of particles, and we call it U(t; 0). 
Then 

T N,t( x >v) ~ <Pt{xWi{v) 

= i<n,Z^(t;0)oJ;o ie W(t;0)n> (1.53) 

+ ^ffi ((O, uHt; 0)aj(«(t; 0) - W(t; 0))O) + (O, («t( t; o) - «t( t . 0))aj,W(t; 0)O) 
y N ^ 

+ ^^((n,Wt(t ; o)a B (W(t;0)-W(i;0))fi) + (O, 0) - Z?(i; 0))a x W(i; 0)O) 

If we define U through 

id t U(t; 0) = (H 2 + H A )U(t\ 0), W(0; 0) = 7, (1.54) 

then it can be proved that 

||W(<; 0)O - U(t; 0)O|| < — = (1.55) 

V AT 



and 

(«(t; 0)Q,MU(t; 0)O) < Ce K *, (1.56) 
which are Lemma 3.9 and Proposition 3.3 of [TU], respectively. Thus, 



Tr 
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V / AT! 4 a /aTT v ' 



Note that the results in [10] are under the assumption that | V(a;) | 2 < D(l — A x ) for some constant 
D. In this paper, we will use these estimates, since the potential V in Theorem 11.11 satisfies the same 
assumption. 

1.5. Factorized initial data. Recall that 

W ] {\fNy)e mNt {a x - y/N<pt(x))e~ mNt W(VNip) = U 1 \t)a x U{i). (1.58) 
By definition, we have that 

( W = a _ N/2 (^/N^f o = 1 (gt(y)) r 

ijv s/Nl 

For factorized initial data, it follows from (ll.58|) and (| 1 . 59[) that 

= (^^aye-^W^fNip)^'^ , (^^a^-^W^fNip)^'^ (1.60) 

'/ v (JV-1) 

(W(v A /vV)W t (t)(a 3 ,+ VNy t {y))U{t)Sl) 

/ , N.(JV-1)\ 

fW'(VA 7 ^)W t (t)(o a! + ViV¥'t(a;))W(t)nj V 
Thus, we obtain the following equation for one-particle marginal. 

d 2 M I /„., rr= „\( Ar - 1 ) f„„ nrz „\( Ar " 1 ) 



N 



^((W(VNcp)li*{t)a y U(t)n) , [W(VN(p)U^(t)a x U(t)nj j (1.61) 



Define 

£(i,z,2/) ~^] t {x;y)-rp t {y) Vt {x). (1.62) 
We also define E2(t,x,y) by putting instead of U in (|1.61l) as follows. 

E 2 {t,x,y) = ^ ^(w(VNtp)ul(t)a y U 2 (t)ny N ~ 1} , (w(VN<p)ul {t)a x U 2 (t)Cl) 

+ ^( 2/ )-^^(PF(v / iV^)fi) : [W(VN^)Ul(t)a x U 2 mj \ (1.63) 

(f2 /, / — , , ^N^- 1 ) /„„, ATT , „\ 
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1.6. Arrangement of the paper. This paper is organized as follows: 

(1) We will first study the evolution of u\{i)a x U 2 {t) in section 2, then, by using it, we will prove 
the following proposition. 

Proposition 1.1. Let E 2 (t,x,y) be defined by (| 1 .63[) . Then, there exist constants C and K 
such that 

Ce Kt 

\\E 2 {t,x,y)\\ Lly < — . (1.64) 

(2) To estimate the difference U{t) — U 2 (t), we study more regularity for U 2 (t) dynamics in section 
3. More precisely, we will prove the following proposition. 

Proposition 1.2. Let ip 2 (t) :=U 2 (t;s)?p be the unique solution of initial value problem 

{idt^Mt)) = N j H 2 {t)ip 2 (t) 

Then there exists a constant C depending only on j, ||Vi||^3, ||V2||l°° and \\<p\\ H i, such that 
\W S Th\\r< C\WH\\ +CT\\{N + 2y +1 ^ 2 \\ F . (1.66) 

(3) Let 

K y (t) :=U ] {t)a y U{t) -U\{t)a y U 2 {t). (1.67) 
In Section 4, using Proposition ll.2[ we prove the following proposition, which gives an estimate 
for the difference between U(t) dynamics and U 2 (t) dynamics. 

Proposition 1.3. LetlZ y {t) be defined as in (|1.67p . Then, for all j G N, there exist constants 
C and K such that 



dy\\{M+l)in v {t)VL\\ 



< 



Ce 



KT 



N 



(1- 



(4) In Section 5, we prove the main theorem, using Proposition ll.il Proposition II .2\ and Propo- 
sition [T31 Some technical estimates are proved in Section 6. 

2. Evolution of u\(t)a x U 2 (f) 

2.1. Properties of U\(t)a x U 2 {t). We first want to obtain some algebraic properties of an operator 
U\(t)a x U 2 (t) , which will act on the vacuum fl. For simplicity, we let 

a x (t) :=u\{t)a x U 2 (t), 

4(«) :=Ut{t)alU 2 (t). 
Formal calculation shows that the time derivatives of a x (t) and a x (t) are given by 



(2.1) 
(2.2) 



Mi 



U 2 (t). 



(2.3) 



Commutators in the right hand side of (|2.3p can be calculated explicitly as follows. 



a x , / dyVya y Vy 



«L / yO.IV ytty 



a x , J dy(V* \<p t \ 2 )(y)ala y ^ = (V * \ip t \ 2 ){x)a x , 
4, f dy{V*\ Vt \ 2 ){y)ala v ] = -(V * \<p t \ 2 )(x)al 
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ax, JJ dydzV ( y ~ z ) L P t (y) L P t ( z ) a ^ a 

4> / / dydzV(y - z)Tp- t {y)ip t {z)ala y 



%>2 /' d y dzV (y ~ z )Vt(y)tpt(z)alal 



Now, (12.31) becomes 



dydzV(y - z)ip t (y)ip t (z)a y a z 



dyV(x - y)Tpi{y)tpt{x)a y , 
dyV(x - y)— t (x)ip t (y)al, 
dyV(x - y)ipt{x)ipt{y)a\, 
dyV(x - y)Tpi(x)Tpi(y)a y . 



a x (t) 



m 

] M ( ~A x a x + (V * \ip t \ 2 )(x)a x 
A x at-(V*\^)(x)at 



U 2 {t) 



U+\( I d V v ( x - y)"Pt(y)^t{x)a v + J dyV(x - y)<p t (x)(p t (y)al 



-Ul{t) 



JdyV(x - y)ip t (x)^t(y)a\ - J dyV(x - y)tp t (x)(p t (y)a y 



U 2 {t) (2.4) 



-A x + (V*\cp t \ 2 )(x) 
dy 







a x (t) 



o A x -(y*|^| 2 )(x) ; V 4(*) 

V(x - y)7pl(y)<pt(x) V(x - y)tpt(x)<p t (y) \ ( a y (t) 



-V(x-y)cpt(x)(p t (y) -V(x-y)<pt(x)ip t (y) J \ a\{t) 
Using (|2.4p . we can immediately see that 



a x (t) 

4(*) 







G 2 \ 






y G3 


G 4 J 


4 



Here, G\, G2, G3, and G4 depend on t, but we omitted it. Since 
we can easily see that G4 = Gi and G3 = G 2 - Thus, 



«(*) 



Gi G2 
G2 Gi 



(2.5) 
(2.6) 
(2.7) 



From (|2.7I) . we can obtain the following lemma, which shows an important property of the operator 
Ul{t)a x U 2 {t). 

Lemma 2.1. Let 

G 2 (t,x) :=Ui{t)a x U 2 {t)VL. (2.8) 
Then, G 2 (t, x) € L 2 (R 3 , L 2 (R 3 )) for any fixed t > 0. Furthermore, if we let 

G 2 (t,x,z) := {Ul(t)a x U 2 (t)Q)(z), (2.9) 
then we have the following bound for G 2 : 

\\G 2 (t,x,z)\\ Llz <Ce Kt . (2.10) 

Proof. We first observe that 

dx(a x (t)n,a x (t)n) = j ' dx(G 2 (t,x),G 2 (t,x)). (2.11) 
On the other hand, we have from the definition and from Lemma 16.21 that 

/ dx(a x {t)n,a x (t)n) = [ dx{a x U 2 {t)tt,a x U 2 {t)n) = {U 2 {t)n,NU 2 {t)n) < Ce Kt . (2.12) 
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Thus, from (f2~TTj) and ([2T2]> . we can find that G 2 (t, x) € L 2 (R 3 , Fx) for any fixed f > 0. In particular, 
Ul(t)a x U 2 (t)tt is well-defined for almost every x. Let r/(") be a function in L 2 (R 3 ") with n > 0. If 
n ^ 1, then the formal calculation (12.7[) shows that 

(^ n \ul(t)a x U 2 (t)n) =0. (2.13) 

Thus, we can see that G 2 (t,x) G L 2 (R 3 , L 2 (R 3 )), which proves the first part of the lemma. Now, we 
rewrite (l2"TTj) and ([2"T2"j) as 

dx(a x (t)Q,a x (t)Q) = J J dxdz\G 2 (t,x, z)\ 2 < Ce Kt . (2.14) 

This completes the proof of the lemma. 

2.2. Proof of Proposition flTTl Using the properties of a x (t), we can prove Proposition ll.il 
Proof of Proposition \l.ll Recall that 

E 2 (t,x,y) 

dl It,.,, r= /„., ATT „\ ( N ~l) 



^nW(VN<p)uUt)a y U 2 (t)nj ,[W(VN(p)Ul(t)a x U 2 (t)nj ) (2.15) 

rfllt , \ (iV-l) / \ (N-l) 

+ Tp t {y)-g= ( (W(VNtp)Q.) , (w{yfN<p)ul{t)a x U 2 {t)SV 



N 

dl //„ w rr= \ ( N ~l) / , 



+ , 

We first want to compute I W(VNcp)K 2 (t)a y Li 2 (t)fl J explicitly. We have already seen from 

Lemma IO that a x (i)H e L 2 (R 3 , L 2 (R 3 )) and we write it as 

a x (t)n=U 2 t (t)a x U 2 (t)n = G 2 (t,x). (2.16) 
Now, we calculate from the definition of the Weyl operator that 

( W ( VNtppl (t )a x U 2 {t) nj 

= e~T (exp(a t (ViV^))cxp(-a(ViV( y 3))G2(<,a;)) (Ar ^ 

= e-*(exp(a t (\/iV^))(G 2 (t,a;)-\/iV : y ( izG , 2(t,x,z)^(z))) (Ar ^ (2.17) 
^^(^(ViV^^^^x) 



TT^yy JdzG 2 (t,x,z)7p(z){aHVNcp)) N l Sl), 



where we used the notation G 2 (t, x, z), which was defined in (|2.9[) . Thus 

(JV-l) 



(W(VN<p)ul{t)a x U 2 {t)n) {x 1 ,x 2 ,--- ,x N -i) (2.18) 

AT-l N-l „ AT-l 

= 3 — 7Tt\ G 2(t,x,Xi) Y[ vixj) — N dz G 2 (t,x,z)lp(z) p( 3 



AT-l N-l n N-l 

We also know that 

AT-l 



(w(VF^)n) (N 1 ( Xu ...,x N - 1 ) = ±l[tp(x j ). (2.i9) 
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Using (|2.17j) and (|2.19j) . we can explicitly calculate the terms in the right hand side of ()2.15|) . The 

first term becomes 

rp I , \(n-i) / \(N-\y 

x v(VNtp)Ul(t)ayU 2 (t)n) ,\W(VNip)Ul(t)a x U 2 (t)n 

N-l N-l „ N-l 



/ (w{\fNip)Ut {t)a y U 2 (t)fi) , (w{VN<pp^t)aMa(t)S^ 

N-l N-l N-l 

= ( E G 2(t,y,x,) H <Pfa)-N J dz G 2 (t,y,z)Tp(z) J[ pfo), (2.20) 

i=l j^i " j=l 

N-l N-l . N-l 

G 2 (t, x, Xi) Y[ <p(xj) -N dz G 2 (t, x, z)Tp{z) Y[ <f{a 



i=l jjii " j=l 

N-l N-l „ N-l 

[X,j 

i=l j^i j = l 

N-l f N-2 f f _ 

= N2 / dz G 2 (t,y,z)G 2 (t,x,z) / dz G 2 {t,y, z)ip(z) / dz G 2 (t,x, z)tp(z). 

Here {■) l 2 denotes the inner product with respect to the variables xi, x 2 , • • • , x^ only. For the second 

X N 

term, we have 

dl //„., r- ,„\( Ar - 1 ) f„ r , f~rz A „\( JV " 1 )' 



^(y)^= UW(VN<p)fl) , {W{VN^)ul{t)a x U 2 {t)Q 

'N-l N-l N-l „ 



3=1 »'=1 ' 2=1 

(iV - 1) y dz G 2 (t, x, z)^(z) -N J dz G 2 (t, x, z)7p(z^j (2.21) 
dz G2(i, a:, z)^(z). 



N 
TV 

= gt(y) 

TV 

A similar calculation shows for the third term that 

Sir, s (N-i) r \ (N-i) 

n(x)^= nw(VN<p)ul(t)a v u 2 (t)n) , [w(Vn<p)si) 

= -^T JdzG- 2 (t,y,zMz). (2.22) 
Therefore, the proposition follows from (f2~2"0"|) , (|2~2T1) . ([2~2"2"|l . and Lemma O 

3. Strichartz Type Estimates for U 2 (t) 

In this section, we prove Proposition 11.21 using Strichartz estimates. For the standard Strichartz 
estimate used in this section, see Theorem 2.3.3 of pQ. 

3.1. Regularity from U 2 (t) dynamics. We will use the notation 

H := J dxV x a) x V x a x (3.1) 

to denote the kinetic energy. 

Lemma 3.1. Let Hq be defined in (|3.ip . Define an operator T on Tx by 

T{t) :=e- ltHo . (3.2) 
Then, there exists a constant C such that, for all tj) G Tx , 

||rWII|<G||^|| 2 . (3.3) 
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Proof. It is trivial that 



|| e -«Ho^(n) (Xn) ||2 =U (n) {Xn)ll 2 (3.4) 



For any fixed i and j, 1 < i < j ' < n, let rj = Xi — Xj. Applying the standard Strichartz estimate for 
Schrodinger operator on we obtain 



\e- itHo ^ n \X n )f 



|| e - 2 «A^(n) (x?i) ||2 



(2f6f2 r 2 



(3.5) 



< C||V (n) PQ|' 2 



17 x^+Xj \{xj ,Xj } X- 



Summing (|3.4I) and (|3.5I) over n, we get the desired lemma. 

Lemma 3.2. Let Hq be defined in ()3.1li . Define an operator $ on JF through 

^{t,X n )= f dse-^-^^^s.Xn) 



for any if> € T . Then, there exists a constant C such that, for all tp G T , 

\\^\\f<CT\\^\\ r , 
Proof. To prove the desired lemma, it suffices to show that 

\\^;\\F n <cT\\^ n \\ T . 

Recall that we only consider a fixed time interval [0,T]. From the standard Strichartz estimate, we 
have 



(3.6) 
(3.7) 
(3.8) 



dse-^-^ H °^ n \ S ,X n ) o <C\\^ n \t,X n )\\ Li 



(3.9) 



For any g(t,X n ) € L^LJ._ X . L Xi+x . \ r x .y we obtain as the dual inequality of Lemma 13. II that 

rT 



dte ttH °g(t,X n ) 







< C\\g(t,X n )\\ L2 e/5 ; 2 



From the standard Strichartz estimate, we also have 

i-T 



ds e isH °^ n \s,X n ) 







<cu^(t,x n )\\ 



(3.10) 



(3.11) 



Combining (|3.10j) and p. lip , we get 



dse~ l{t - s)Ha ^ {n) {s,X n ) 



sup / dt j ds(g(t,X n ),e-^-^ H °^ n Hs,X n )) 



(3.12) 



9 Jo Jo 

rT 



supM dte ttH °g(t,X n ), j dse isHo i>W( S) X n ) 







<Csup|| 5 (t,X n )|| L2L a/5 L2 L2 =C||V (n) (*,X n )|| i i J 



where the supremum is taken over all funtions g satisfying 



\\g(t,X n )\\ L2L e/ 5 



(3.13) 
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Applying Christ-Kiselev lemma to (|3.12j) . we finally obtain that 



ds e- l(t - s)Ho ^ n) (s,X n ) 



T 2 T G Tl T2 



<c\\^ n \t,x n )\\ Ll 



Thus, it follows from (|3J]) and (|3~T4l) that 



||*i n) ||F„ < C\\^ n \t,X n )\\ LiL , < CTU^(t,X n )[ 



(3.14) 



(3.15) 



which proves the claim (|3.8|) . This concludes the proof of the desired lemma. 

3.2. Estimate on H^it)- To prove Proposition [L^J we need the following lemma. 

Lemma 3.3. Let V = V\ + V% with V\ € L 3 and V2 S L°° . Then, for all j g N, there exists a constant 
C depending only on \\Vi\\l3, ||V2||z,°° and \\ip\\iji, such that 



\\Af j (H 2 (t) ~ Bo)1>\\ < Cll + 2) 3+ Vll 



for all if) G Tx ■ 
Proof. Let 



L := dx{V * \(p t \ 2 )(x)ala x , 



M 



B 



dxdyV(x - y)if t (x)ip t (y)ala x , 
dxdyV(x - y)Tp- t {x)Tpl{y)a 



Then, H 2 - H = L + M + B + Bi. 
Note that 

in * M 2 L < U^MMUM^ + Mi-ll* < a 

Since L commutes with Af and V * \ft\ 2 € L°°, there exists a constant C such that 

\\Af 3 L^\\ = \\LAf^\\ < CH^'+Vll- 
To prove an estimate for ||A/' : 'M-0|| = jjM.APV'll ? we observe that, for any £ e Fx-, 



\(£,MAf^)\ 



dxdyV(x-y)ip t {x) V t{y){{N +l)- 1/2 a y ^{N +l) 1/2 a x N^) 



dy\\{M+l)- l > 2 a y £,\\ 



1/2 



1/2 



(3.16) 



(3.17) 
(3.18) 
(3.19) 

(3.20) 
(3.21) 



(3.22) 



<cneiiii^' + Vii- 

Since £ was arbitrary, this shows that 

\\M j Mip\\ < C\\M j+1 i)\\. (3.23) 
Estimates for \\Af j Bip\\ and HAP'-B^II 

can be obtained similarly. Since 

dxdy \V(x - !/)| a Ma:)| 2 M!,)| 2 < IN|| 3 |M|le + H^Ul- |M|£ 2 |M|1, 2 < C, (3.24) 
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we have that, for any £ £ Fx, 



\{Z,WB4)\ 



dxdy V(x - y)(p t (x)<p t {y)(€,a x a y (Af - 2) 3 <ip) 



1/2 / , \ 1/2 

< ( / / dxdtf -»)| a |^(a:)| 2 |^(y)| 2 ||e|| 2 ) / dxdy\\a x a y (M-2y^\\ 2 



< C||£||||AP' + V||- (3.25) 
Again, since £ was arbitrary, this shows that 

\\M j Bil)\\ < C\\Af j+1 i()\\. (3-26) 

Similarly, we also have that 

\\Af j B^\\ < C\\(Af+2y +1 i/j\\. (3.27) 
Thus, from (|3~2T|) . (|3~23| . I|52B|) . and (pH?) . we get 

||AP'(tf 2 (t)-tfo)?A|| < \\AT ] L^\\ + \\AT 3 MtP\\ + WW (B + B^)i;\\ < C\\(N + 2) j+1 ipl (3.28) 
which proves the desired lemma. 

3.3. Proof of Proposition 11.21 From Lemma l3~Tl Lemma I3T21 and Lemma l3~3l we can prove Propo- 
sition 11.21 

Proof of Provosition \1.SX By using the group e~ ltH ° generated by Hq, we can write the evolution 
of7V J > 2 (t) as 

M j ^ 2 {t) = e- ltHo Af j ip+ f dae- t{t ~^ Ha M j {H 2 {o-)- H )ip 2 {o-). (3.29) 



Thus, by Lemma 13. 11 Lemma l3.2[ and Lemma 13.31 we have 

WM?< C\\M^\\+CT\\^{H 2 (-)- Ho)M-)\\? (3-30) 

<c\\Wip\\+CT\\(Af + 2y +1 ^ 2 y, 

which was to be proved. 

4. Difference between U(t) Dynamics and U 2 (t) Dynamics 

4.1. Estimates on Hs(t) and H4. From the regularity we have seen in Proposition ll.21 we can obtain 
the following estimates. 

Lemma 4.1. Let V — V\ + V 2 with V\ E L 3 and V 2 £ L°° . Then, for all j G N, there exists a constant 
C depending only on j , ||Vi||l3, ||V2||l~ and \\(p\\hi, such that 

WH z {m\<-%w+iy +i n (4.1) 

v jy 



for all tp £ Fx- 
Proof. Let 

M{t) = JJ dxdyV(x - y)Tp t (y)ala y a x . (4.2) 

Then, 

Af>H 3 (t) = -L(AP"4(i) +WA 3 (t)). (4.3) 
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Now we estimate Afi A' 3 (t) and Af^A 3 {t) separately. The first term Af^A\(t) satisfies 

\w^A\m{x n )t L ^ 

= ^^' +1 (n-l) 2 / dx 3 ---dx n dx 1 dx 2 \V(x 1 - x 2 )\ 2 \ l p t (x 2 )\ 2 \^ n - 1 \x 1 ,x 3 , - ■ ■ ,x n )\ 2 

n=3 J J 

+ V+ 1 J dxidx 2 \V{x l -x 2 )\ 2 \y t {x 2 )\ 2 \il>V{xi)\ 2 (4.4) 

OO 

<2^^ +1 (n-f) 2 



n— 3 



>%n)\\h 



+ 2^ 2 (||^|||3||^||| s + ||F 2 |||^ t ||^IIV> (1) (*l)||i= 

x 1 

<C||(AA+l)^+i^|| 2 , 
and Afi A 3 (t) satisfies 
\\^A 3 (t)^(X n )f L2xn 

= jr n i+ 2 (n + l) fd Xl ---dx n I dyV{x 1 -y)ip t {y)^ n+1 \y,x 1 ,--- ,x n ) % (4.5) 
n=l J ^ 

<2^n^ 2 (n+l) / ^•••da ;n (||F 1 ||| 3 ||^||| 6 + ||y 2 ||l oo ||^||i0ll^ (n+1) (y^i,--- 
n=l ^ 

Hence, from (|4~3]) . (|4~4|) . and ({475]) we get 

\\Af j H 3 (t)^\\ < -^=\\{M+iy + Hl (4.6) 
V TV 

which was to be proved. 

Lemma 4.2. Let V = V\ + V 2 with V\ £ L 3 and V 2 £ L°° . Then, for all j £ N, there exists a constant 
C depending only on j , || Vx || 2^3 , ||V^||j,oo and ||^|| jyi , such that 



C CT 

d s \\Ar^H 4 U 2 ( S ; Sl )^\\ 2 < -^||^+ 2 ^|| 2 + ^\\(Af+2y+ 3 U 2 (.-,s 1 )i;\\ 2 r (4.7) 



for all ij) £ Tx ■ 
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Proof. From Holder's inequality, we have 

n J ; ^ A 



n=2 



*<3 



< — V n 2j 
N 2 2—i 

n=2 



n(n — 1) 



n . 

i<j ■> 



n=2 i<j J 



x / + Xj)dxi ■ ■ ■ d,Xi ■ ■ ■ dxj ■ ■ ■ dx n (U2(s;si)ip) (xi,---,x n ) 



- ■— > i<j n 
+ \\V 2 \\l^\\{U 2 {s ]Sl )^) {n) (x ir -- ,s„)lli»J- 



n=2 



C CT 

^ ds\\WH 4 U 2 (s; Sl )yj\\ 2 < ^||^ +2 W 2 (-; S i)VII| + ^ll^ +2 W 2 (-; S i)VII^ 



Thus, 



Proposition II .21 shows that 

||AP"+ 2 W 2 (.; Sl )V|| j? < C||^'+V|| + CT|| (TV + 2)->+ 2 W 2 (-; 

Therefore, 



f 1 r*T 
d.s||^iJ 4 W 2 ( S ; Sl )^|| 2 < _||^+V|| 2 + — m + 2y+ 3 U 2 (- :Sl )^, 



N 2 



which was to be proved. 



(4.8) 



(4.9) 



(4.10) 



(4.11) 



4.2. Proof of Proposition [TT3l We are ready to estimate the difference between W (t)a y U(t)Cl and 
U\{t)a y U 2 (t)Vl by proving Proposition 11.31 Note that we are free to use the Duhamel formula with 
the operators U and U 2 on a dense subset V(Hq) n T>(M) of Fx, which can be seen from equations 
fL"4"5|l and (|T3T|) together with the fact that the Weyl operator W(y/N(p a ) maps V(H ) onto V(H ) 
if (p s € H^R 3 ). (See also Lemma 3.1 of @].) 
Proof. Let 



and 



Then, n y (t)=nl(t)+TZ 2 y (t). 
Since 



H\{t) := (U\t)-Ul(t))a y U(t) 



H\{t) :=Ul{t)a y (U(t)-U 2 {t)). 



(4.12) 



(4.13) 



U\t)-Ul(t)=i / dsU\s;0)(H 3 {s)+H 4 )ul(t;s), 



(4.14) 
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applying lZ y (t) on the vacuum gives 

+ l)in l y {t)Q\\ = \\{M+l)i{U\t) -Ul(t))a y U(t)n\\ 

(4.15) 



ds{N + l)4w t (a)(£T 3 (s) + H 4 )U 2 (s;t)a y U(t)n 

< I ds\\(Af + l)itf(s)(H 3 {s)+H4,)U2(s;t)a v U{t)n\\ 
Jo 

= [ ds(U\s)(H 3 (s) + H i )U 2 {s;t)a y U(t)n,(Af+iyU f (s)(H 3 (s) + H 4 )U 2 (s;t)a y U(t)n)? . 
Jo 

Thus, from Lemma 16 . 1 1 with Schwarz inequality, we obtain that 

dy\\(Af+ l^TZlmW 2 



<Ct I dy I ds\\(Af + l) j+1 H 3 (s)U 2 (s;t)a y U(t)Cl\\ 



(4.16) 



Now, from Lemma 16.11 Lemma 16.21 and Lemma 14.11 the first term in the right hand side of (|4.16p 
can be estimated as 



+ Ct dy ds\\(Af + l) j+1 H 4 U 2 (s;t)a y U(t)n\\ 



dy / ds\\(Af +iy +1 H 3 (s)U 2 (s;t)a y U(t)n\\ : 



< 



C_ 
N 



dy / ds\\(Af + iy + ^U 2 (s;t)a y U(t)n\\ : 



< 



< 



Ce KT 

N 
Ce KT 



dy \ ds(a y U(t)n,(AT + l) 2j+5 a y U{t)n) 



(4.17) 



< 



N 
le 1 



ds(U(t)tt, (TV + l) 2j+6 U(t)n) < 



Ce KT rt 



o 



N 



ds{fl, (Af + 1) 4 ^+ 14 0) 



Ce KT 



N 



Similarly, from Lemma 16.11 Lemma 16.21 and Lemma 14.21 the second term in the right hand side of 
(|4.16p can be estimated as 



dy / ds\\(Af+ iy +1 H 4 U 2 (s;t)a y U(t)n\\ 



< 



C_ 



dy\\(Af+iy +3 a y U{t)n\\ 



CT 



dy\\(Af +iy +4 U 2 (-;t)a y U(t)n\\ 2 jr 



r Cp kt 

n p KT Cp kt 
Hence, from (14TTH)) . (|4~T71) . and (14TT51) . we get 



(4.18) 



d ! /||(.A/'+l)*ftJ(t)fi|| 



< 



Ce 



KT 



N 



(4.19) 
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The study of 1Z 2 (t) is similar and gives 



dy\\(Af + l)ill 2 y (t)nf < 



Ce 



KT 



Therefore, 



which was to be proved. 



|(7V+l)»7^(t)Q|| 



< 



N 

Ce KT 



N 



(4.20) 



(4.21) 



5. Proof of Main Theorem 



Now, we are ready to prove the main theorem. 

Proof of Main Theorem. First, we write the difference between E(t, x, y) and E2(t, x, y) in detail as 
follows. 



E(t,x,y) - E 2 (t,x,y) 



r 2 

N 



(N-l) 



" :x / (w(VN<p)n y (t)n 

(w{\fNip)Ul{t)a v U 2 {t)tt 

(w(VNip)n v {t)Q. 
d 2 

N 
d 2 



N 
N 



, [W{\/N(p)ul{t)a x U-2(t)n 

(N-l) 



(w(VN(p)n x {t)n 
Wiv)^= ^(w(VNv)n) ( "~ L \ (w(VN<p)n x (t)n 

d 2 / / , \ (N-l) / , 

+ ipt(x)^nw(VN<p)K y {t)n) ,{W(VN(p)fl 



(N-l) 
(N-l) 



(5.1) 



+ 



(N-l) 
(N-l) 



(N-l) 
(N-l) 
(N-l) 



'N 

From Lemma 16.21 and Proposition 1 1 . 31 the first term in the right hand side of (|5.1[) can be estimate 



as 



f I / </■<■".</ 



(N-l) 



< 



< 



U N 

N 



l N 



dy\\K y (t)n\\ 



dx\\l4(t)a x U 2 (t)n\\ 



(5.2) 



Ce 



KT 



Ce KT < 



Ce 



KT 



Similarly, the second term can be estimated as 



N 



dxdy 



{w{VN v )ul{t)a y U 2 {t)d s j {N 1] , (w(y/N(p)Hx(t)!l 



(N-l) 



2\ 2 Cp KT 

<%-■ (5-3) 



To estimate the third term, we use the following estimate 

(N-l) 



N 



dxdy 



< -f ( / dy\\n y m\ 



(w(VN(p)n y (t)S}) ,(w(VNip)Tz x (t)n 



(N-l) 



2\ 2 



dx\\n x {t)n\\ 



< 



NvN 



(5.4) 
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To estimate the fourth term in the right hand side of (|5.1j) . we hrst note that 



(N-l) 



= /^ (Ar - 1} , (w(VNip)n x (t)n 



(N-l) 



^(ip 9 &- 1 \W(VN<p)Tl x (t)n) = ^(w(v / iV^)V 8(Ar " 1) ,^x(t)0) (5.5) 



It follows from Proposition 11.31 and Lemma 15731 that 



dx 



< 



< 



Ce KT 



dx 



(AT + l)*1l x (t)Sl 



~ d N VN 

Thus, from (|5 .5[) and (|5.6p . we obtain 



^xdy ^(y)-^/(w(v^V^)0) (iV X \ (w(VN<p)K x (t)fL 



(N-l) 



< 



Ce 



KT 



N 



The last term in the right hand side of (|5.ip can also be similarly estimated as 



< 



dxdy 

Ce KT 
N ' 



(N-l) 



(N-l) 



2\ 2 



(5.6) 



(5.7) 



(5.8) 



Therefore, together with Proposition II .11 inserting (|5.2I) . 
yields 

\\E(t,x,y)\\ Lly <\\E 2 (t,x,y)\\ Lly + 



153), and, g3) into (ED) 

(5.9) 



< 



7Y ~ 7Y 

By Remark [2] and Remark [3j this completes the proof of the main theorem. 

6. Expectation of Number Operator with respect to the Various Evolutions 

The following lemma shows that the expectation of the number operator with respect to the full 
evolution is bounded uniformly in N. 

Lemma 6.1 (Proposition 3.3 of [ID]). Suppose that V = Vi +V 2 where V% £ L 3 (M. 3 ) and V 2 £ L°°(R 3 ). 
LetU(t;s) be the operator satisfying (11.43[) . Then, for all j £ N, there exist constants C and K, 
depending only on j , \\ip\\ H i, and ||V||l 3 such that 

(U{t;s)^,AfiU(t; S )i>) < Ce K ^ (V, (/V + l) 2j+2 ^), (6.1) 

for all ip £ Tx- 
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Remark 5. Specially, ifip = Q., we have 

(U(t; s)fi, M 3 U(t; s)Q) < Ce Klt - sl . (6.2) 

A similar result holds for lA 2 (t) evolution. For completeness, we prove it here. 

Lemma 6.2. Suppose that V = Vi + V 2 where Vi G L 3 (R 3 ) and V 2 G L°°(R 3 ). Let U 2 (t; s) be the 
operator satisfying (jl.Sip . Then, for all j G N 7 there exist constants C and K, depending only on j, 
and ||V||£ 3 such that 

(U 2 (t;s)xj;,Af j U 2 (t;s)i(}) < Ce Klt ~ sl (V, (Af + 1)», (6.3) 

for all %l> G Tx- 

Proof. It suffices to prove the lemma when s = 0. We have 
d 



dt 



(U 2 (t; 0)^,(Af + l) J U 2 (t;0)^) = (U 2 (t;0)ip, [iH 2 , (Af + iy]U 2 (t;0)ip) 



= Im ff dxdyV(x~y)Mx)My)(U2(t;0)4>,latal(Ar+iy]U 2 (t;0)iP) 



3-1 



k=0 



dxdyV(x - y)<pt(x)ip t (y) 



(6.4) 



x ((a x (Af + l)*M 2 (t; 0)V», a\(Af + 3)*H 3 (t; 0)V>) + (a x Af?U 2 (t; Q)i/>, a\(Af + 2)§U 2 (t; 0)V>) 
j'-i 



fc=0 



sup||F(x--M-)ll 



L' 2 



dx\ Vt (x)\\\a x (N + l)*U 2 {t; 0)^|| || (Af + i)—U 2 {t; 0)V|| 



dx\ipt(x)\\\a x N 2 U 2 (t; Q)tp\\\\(Af + 3)~U 2 (t; 0)ip\\ 

< C\\ (M + l)iU 2 (t; 0)V>|| 2 = C(U 2 (t; Oty, (A/ - + l) j U 2 (t; 0)^). 
Since ^(0; 0) = /, we also have 

(W 2 (0;0)V,(AT + 1)^ 2 (0;0)V) = (^,(AA+1)». 
Using (|6.4|) and (|6.5[) the conclusion follows directly from the Gronwall's lemma. 
Lemma 6.3. There exists a constant C such that 



{AT+ l)-^W{^fN v )\v® {N ~ 1] ) 



< 



C 



<N 



(6.5) 



(6.6) 



Proof. It is obvious that W(v / AV) t ((p^ N ~^) is a linear combination of tensor products of ip. We 
define the coefficients A m such that 

(m) 



For m < N — 1, we can explicitly calculate A m - Since 



(^(Va^^- 1 ))) ™ = ( e - JV / 2 exp(a t (-ViV^))exp(a(v / 7V^))(^( A '- 1 ))) 



we have 



m I _ 

^rj(-v / iV)' t Vm - + lVm - k + 2 • • 



(6.7) 

(6.8) 
(6.9) 



(A - m - 1 + fc)! 



■N-m—l+k , , , 

A v / A-l v / A-2 ■ ■ • y/m-k + l. 
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We want to estimate \A m \. When m = 0, we have 



N-l 



An = er N/2 - , 

y/(N-l)\ d N 

When 1 < m < N — 1, A m can be simplified as follows. 



A m = e-^VN N - m - 1 J^^J2(-l)^ 



-N/2 



N 



N—m—l 



-L 



(JV-m-l) 



(m - k)W.(N - m - 1 + k)\ 



(6.10) 



(6.11) 



V (iV-l)! 

where the associated Laguerre polynomial L„ (x) is defined by 



L^(x) ■= Vf-l) fe (" + ") 



k=0 



k\(n - k)\{a + k)\ 



(6.12) 



We now follow the argument used in the proof of Lemma 4.2 in [10]. The sharp estimate obtained 
by Krasikov [5] for iffi (x) shows that 



\LM(x)\ < 



(n + a)\ x(s 2 — q 2 ) M c+ 



r(x) 



e 2 x 2 



(6.13) 



where 

s = (n + a + 1) 2 + n» , q = (n + a + 1)3 - , r(x) = (x - g 2 )(s 2 - x). (6.14) 
Thus, we obtain that, for 1 < m < N — 1, 



| An | < e-^W 

< civ- 1 /^- 1 / 4 . 



(iV-1)! 



JV-m-l / m! (N-i)\ / 4^^^ 



i! V 47Vm - 



3 iV/ 2Ar 



4%/iVr 



4Nm — in 



For m > N, we only need the following bound. 

oo 

E \A m \ 2 <\\(w(VN^(^ N -V)) 

m=N 

Therefore, from f)6. 10[) . (|6.15[) . and (|6 . 16[) . we obtain 

(^+i)-^(v / iv^ t (^ (jv - i) ) - y 

z — ' 771+1 



< 1. 



(6.15) 
(6.16) 

(6.17) 



JV-1 



C 



CN ~ 112 ^ J^TW' + ^ z ~^ l " 4 '"' 2 - CN ~" 2 - 

Jv m— 1 v ' m—N 



This proves the desired lemma. 



7. Stability of the Operator 



The following lemma shows that the operator H2(t) is stable. (See Proposition 3.4 of [7] for more 
details.) 

Lemma 7.1. There exists a constant C, K > such that, for the operator ^(t) = H2(t) +C(jV+ 1), 
we /lave i/ie operator inequality A%(t) < KA2{t), where A2{t) = (d/dt)A2(t). 
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Proof. Let 

dH 2 (t) 



H 2 (t) := 



dt 



= JJ ' dxdyV(x-y)<pt{y)<Pt{y)aia x + J J dxdyV (x - y)ip t {x)w t {y)a\a x 

+ JJdxdyV{x-y)ipt{x)ip t {y)a\.a\ + h.c. (7.1) 

where h.c. denotes the Hermitian conjugate and ip t = d t (p t . 

In order to control (tp, H 2 (t)ip) for e Tx, we need to estimate terms such as 



dxdy{ip, V(x - y)(pt{x)ip t {y)a\alip) . (7.2) 
We know that ip t is the solution of the nonlinear Hartree equation and satisfies 

My) = -i[-*Mv) + (v * \<p t \ 2 ){v)Mv)]- (7-3) 

Thus, 

J J dxdy(i>,V(x-y)<pt(x)<pt(y)alali>) 

dxdy (a x a v ^ 7 V(x - y)y t {x)[{- Ay t ){y) + {V * \<ptf)(v)Mv)]il>) ■ ( 7A ) 



Since (p t € H 1 and 

\\V*\<p t \ 2 \\co<C, (7.5) 

we obtain that 



J J dxdy(a x a y ip,V(x - y)tp t (x)(V * \ip t \ 2 )(y)(p t (y)ip) 



OO 

< £ V^^T)H (n - 2) \\ ■ U {n) \\ ■ \\v(x - y)vt{x){v * \<Pt\ 2 )(y)vt(y)hi y (7.6) 



n=2 



By integrating by parts, we find that 

J J dxdy (a;, y, X n _ 2 )V{x - y)<p t (x)(—Aip t )(y) 

= JJ dxdy ({y v ^)){x,y,X n _ 2 )V{x - y) - ^)(x,y,X n _ 2 )(VV)(x - y)) • <pt(x)VMv) 
= JJ dxdy (V y ^ n ))(x,y,X n -2) ■ V(x - y)<p t (x) V<p t (y) (7.7) 
+ JJ dxdy (-V x ^)(x,y,X n _ 2 ) ■ V(x - y)y t (x)V <p t {y) 
f I dxdy ^ n ){x,y,X n - 2 )V{x - y)(-V<pt)(x) • Vp*(y). 
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Thus, we get 



dxdy (a x a y ip, V(x - y)(p f (x)(-Aip t )(y)ip) 



(7.8) 



= | Vn(n-1) / dxdydX n _ 2 ^){x,y,X n _ 2 )V{x - y)^ t {x){-^ t ){y)^ n ' 2) {X n ^ 2 ) 

n=2 
oo 

< J2 - 1) 



n=2 
X 



(\\(V v ^)(x,y,X n _ 2 )\\ L , \\v(x- y ) Vt (x)vMy)\\Li v H {n - 2) (x n _ 2 )\\ L2 

+ \\(\7 x ^)(x 7 y,X n _ 2 )\\ L 2 \\V(x-y)ip t (x)VMy)\\Ll (X n _ 2 )\\ L% 

x,y.X n 2 x ;V X 

+ \\V(x-y)^ n Hx,y,X n ^ 2 )\\ L 2 \\V<Pt(x)VMv)\\ii H (n - 2) (X n _ 2 )\\ L , 

When V 2 < D(l - A), 

\\V{x- y )^ n \x, y ,X n _ 2 )\\ L . 



<D(\\(V x ^)(x,y,X n _ 2 )\\ L 2 +U (n \x,y,X n ^ 2 )\\ L 2 
Hence, from (|7.8p and (|7.9p . we have 

dxdy (a x a y ip, V(x - y)(p t (x)(-A<p t )(y)ip) 

O 

<cV i)||(v^ (n) )(^y,x n _ 2 )|| i2 ||^ n - 2 )(x„_ 2 ) 

y r a;,y.A TI _2 



n=2 

OO 



+ CJ2 Vn(n - l)U (n) (x, y, X n _ 2 )\\ L , \\^ n - 2 \X n _ 2 ) 

n=2 

<C(i>,(H +Af + 
Now, ([Til]) . (EH), and (fTTTOt show that 

dzdtfM V(z - 2/)^(a;)^(y)4 a t^)| < C(i/>, (H + C(Af+l))1>). 

Other terms in the right hand side of (|7.1|) can be estimated similarly. Hence, we get 

Uj,H 2 m) < {H + C(7V + 1))V>. 
Furthermore, we also know that H 2 (t) — Ho > —C(N +1), or 

flaW + CCAT+l) >H Q . 
(See Corollary 2.1 in [5].) Then, we can find constants Co and C\ such that 

j t (H 2 (t) + C (Af + 1)) V) < (ff a (*) + 0>(AA + 1)) V)- 
This proves the desired lemma. 
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